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We study transport properties of a strongly interacting superfiuid mixture of two Fermi-liquids. A 
typical example of such matter is the neutron-proton liquid in the cores of neutron stars. To describe 
the mixture, we employ the Landau theory of Fermi-liquids, generalized to allow for the effects of 
superfluidity. We formulate the kinetic equation and analyze linear response of the system to vector 
(e.g., electromagnetic) perturbation. In particular, we calculate the transverse and longitudinal 
polarization functions for both liquid components. We demonstrate, that they can be expressed 
through the Landau parameters of the mixture and polarization functions of noninteracting matter 
(when the Landau quasiparticle interaction is neglected). Our results can be used, e.g., for studies 
of the kinetic coefficients and low-frequency long-wavelength collective modes in superfiuid Fermi- 
' mixtures. 
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I. INTRODUCTION 



Neutron stars are ultracompact objects with radius R ~ 10 km and a mass ~ M Q , where M Q « 1.99 x 10 33 g is the 
solar mass. The density in the cores of neutron stars is a few times larger than the density po ss 2.8 x 10 14 g cm -3 in 
atomic nuclei. The simplest matter composition in the internal layers of neutron stars includes neutrons, protons, and 
electrons. Neutrons and protons form a strongly nonideal, degenerate Fermi-liquid, while electrons can be considered 
■ as a degenerate gas of free relativistic particles. It is generally accepted [l|-[3| that at low enough temperatures 
^p H ' T < (10 8 — 10 10 ) K neutrons and protons become superfiuid. Thus, to compare the theory with observations of 
neutron stars one should be able to describe a strongly nonideal, superfiuid mixture of Fermi-liquids. In particular, 
' the kinetic properties of such matter are of primary interest. 

One of the important problems in the kinetics of superfiuid mixtures is the calculation of the so called response 
functions that characterize response of the system to small external perturbations. The response functions determine, 
for example, the emissivity of various neutrino processes (e.g., the Cooper- pairing neutrino emission processes 04^1 or 
t-H ■ bremsstrahlung processes [10|]) and thus influence the thermal evolution of neutron stars. Furthermore, the formalism 
of response functions is ideally suited to study collective modes in the system which may lead to a discovery of new 
exotic processes of neutrino emission. Finally, the vector response functions determine the screening properties of 
*~Q ' particle interaction, influence particle collision amplitudes and hence affect kinetic coefficients. 

Although for normal matter of neutron stars the response functions and the collective degrees of freedom have been 
extensively discussed in the literature (see, e.g., Refs. [Tl| - |22l |). the author is unaware of any discussion concerning 
the response in superfiuid strongly interacting mixtures. However, for a one-component strongly nonideal superfiuid 
Fermi- liquid (neutrons in the crust of a neutron star), the response functions were calculated recently by Leinson [9|. 
He generalized the results of Leggett [23J to the case of not too small wave vectors and perturbation frequencies and 
calculated the response of the system to axial-vector perturbation. 

The same problem was also examined in Ref. @. The authors of this reference considered a superfiuid neutron- 
proton mixture; however, they do not allow for interaction between neutrons and protons. Therefore, de facto they 
analyze the system of two noninteracting one-component Fermi-liquids. The results obtained in Ref. [(| were discussed 
(and criticized) in Ref. 

It is convenient to calculate the response functions using the method of kinetic equation. For that one has to 
formulate the kinetic equation which correctly describes (a) possible superfluidity of neutrons and protons and (b) 
various Fermi-liquid effects connected with interactions between neutrons and protons. For a one-component Fermi- 



liquid such an equation was derived, for example, by Betbeder-Matibet and Nozieres [24| (see also Ref. 1251 1 ) . N eglecting 
the Fermi-liquid effects the superfiuid kinetic equation was thoroughly examined, e.g., in Refs. [26H30l | (see also 
references therein). 

There exist a number of papers exploring the transport properties of superfiuid matter in the cores of neutron stars 
(see, e.g., Refs. [31-35]). However, to our best knowledge, the kinetic equation, satisfying both conditions (a) and 
(b), has not yet been proposed. 

Partially, this is because most studies in the literature have been devoted to kinetics of electron gas, in particular 
to electron kinetic coefficients. To calculate these coefficients the kinetic equations for nucleons are not required since 
one usually neglects the interaction between neutrons and electrons and considers protons only as scatterers. Yet 
another reason is the one-component Fermi-liquid (for example, electrons in metals or liquid helium-3), for which the 
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Fermi-liquid effects are well known to play no major role in a variety of kinetic phenomena. 

Here we demonstrate that the situation with the superfluid mixtures is different and that the Fermi-liquid effects 
can be more pronounced there. To a large extent this is due to the mean-field interaction between neutron and proton 
quasiparticles. In the hydrodynamics of superfluid mixtures [3(l - l4lj this interaction leads to entrainment of neutrons 
by the superfluid motion of protons and vice versa (the so called entrainment effect). We will show that the similar 
effects are important for the kinetics of superfluid mixtures. 

In the present paper we formulate the collisionless kinetic equation, satisfying both conditions (a) and (b). Em- 
ploying this equation, we calculate and analyze response of the system to a vector (electromagnetic) perturbation 
(more precisely, we calculate the longitudinal and transverse polarization functions). To describe the superfluid mix- 
ture we use the framework of the Landau theory of Fermi-liquids, extended by Larkin and Migdal [42l l43j and by 
Leggett [23L |44| to allow for superfluidity. For simplicity, we assume that in thermodynamic equilibrium both particle 
species pair in the spin-singlet state (this is indeed a simplification since, according to microscopic calculations 
|2|, neutrons in the cores of neutron stars pair in the spin-triplet 3 P2 state). 

In the collisionless approximation, which we are interested in here, the electrons interact with the nucleons only 
through a self-consistent electromagnetic field. Thus, in what follows, the electrons can be safely ignored; this does 
not change our results. 

The paper is organized as follows. In Sec. II we derive the kinetic equation describing, in the linear approximation, 
mixtures of superfluid Fermi-liquids. In Sec. Ill we calculate and analyze the longitudinal and transverse polarization 
functions for mixtures. In Sec. IV we propose a nonlinear kinetic equation describing superfluid mixtures in the 
quasiclassical limit. Section V presents the summary. 

Below, unless otherwise stated, we use the system of units in which the Planck constant h, the speed of light c, the 
Boltzmann constant k^ , and the normalization volume V equal unity, H = c = fes = V = 1. 



II. THE KINETIC EQUATION FOR SUPERFLUID FERMI MIXTURES 
A. Thermodynamic equilibrium 

To establish notations, let us briefly consider a strongly interacting degenerate homogeneous Fermi-liquid in thermo- 
dynamic equilibrium at some temperature T . We assume that the liquid is two-component and composed of particles 
of two species, i = 1 and i = 2. Here and below the indices i and j refer to these species. 

Weakly excited states of our system can be described in terms of quasiparticles. According to the Landau theory 
of Fermi-liquids, in the normal (nonsuperfluid) matter, the energy of a quasiparticle i with momentum A; is a 
functional of quasiparticle distribution functions A^^ 1 (see, e.g., Refs. (45l. l46j). 
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For nonsuperfluid matter in thermodynamic equilibrium the distribution function A/^ is denoted n^ and is given 
by 



(2) 



In Eqs. (JTJ and © e^ Q is the quasiparticle energy at T = 0; a and a' = ±1 are spin indices; 0^ = 9(kpi — k), 
where 6(x) is the step function, and k^i is the Fermi momentum. Furthermore, a^ a = ajS or ajj?? is the annihilation 
operator of a quasiparticle in a state (ka) and /i.j is the chemical potential. Finally, / y (fe, k ) in Eq. ([TJ is the Landau 
quasiparticle interaction. In the present paper we only deal with the spin-unpolarized matter and do not consider 
forces that rotate spin. This allows us to disregard the spin dependence of the interaction function, assuming that 
f l j(k, k') is already averaged over spin variables. We also suppress spin indices in formulas, whenever possible. Notice, 
however, that the kinetic equation (|19|) . obtained below, can be used, with minor modifications, to study transport 
properties of spin-polarized matter [then, of course, one should take into account the spin-dependence of f^(k,k')]. 

In the vicinity of the Fermi surface the lengths of the vectors k and k' in the arguments of the function ) 
can be approximately set equal to k ~ k-pi and k' ~ kpj, while the function itself can be expanded into Legendre 
polynomials P; (cos 6) , 

/«(*,*') =j2fi 3 p i ( cos 0)> (3) 
1 
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where 9 is the angle between k and k' , and j\ 3 are the symmetric Landau parameters, = ff 1 , 

Assume now that both particles, i — 1 and i = 2, are superfluid (and we are still in thermodynamic equilibrium at 
a temperature T). The spin 1/2 elementary excitations in superfluid Fermi- liquid are the Bogoliubov excitations. In 
the absence of superfluid currents in the system, their energy £^ and the distribution function Jj*q are (see, e.g., 
Ref. [13) 



where £^ equals (see, e.g., Ref. [38| 
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jUi = v Fi (k - fe Pi ) + O 



(6) 



The function A/^q will be defined in Eq. p^|) ; vpi = fcpi/m* is the Fermi velocity and m* is the effective mass. In the 
nonrelativistic limit m* can be found from the equation [33, [H, [48[ 
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where 

1 " 3tt 2 

is the number density; is the bare mass; and the symmetric matrix Gij is defined as 

G y = T^PfiPfj fi ■ 

Furthermore, Aj in Eqs. (U])-© is the energy gap in thermodynamic equilibrium. Because we assume the singlet- 
state 1 Sq pairing of quasiparticles, the gap Aj depends only on k — \k\. It can be found from the standard equation 



(8) 



(9) 



A i (k) = -Y,V {l) (k,k') Ai(k') F$ 
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Here we define 



F ® = 



2E, 



(i) 



tanh 
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(11) 



In Eq. @0J) V®(k,k') is the pairing potential for particles i. In analogy with the function / SJ (fc,fc ) [see Eq. {3J] 5 we 
expand it into Legendre polynomials, 



V (l) (k,k') = J2 V i (l) F '( cos ' 



(12) 



Near the Fermi surface a smoothly varying function Aj(A;) can be approximated as Aj(A;) w A^Api) (see, e.g., Ref. 
'24]). Then, combining Eqs. (|TU|) and ((HJ), one obtains the following equation for Aj, 



i = -k.™ E p y 



(0 



(13) 



As follows from Eq. ([S]), the distribution function for Bogoliubov excitations is a scalar quantity. In contrast, 

(i) 

the distribution function nj.^ for Landau quasiparticles in superfluid matter is a matrix rather than a scalar. In 
thermodynamic equilibrium it can be written as (see, e.g., Ref. p4| ) 
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where i is the unit matrix. The 'energy matrix' ej^ in Eq. (| 14[) equals 

(15) 



c fco-0 



<4 



r(i) 

For superfluid matter in the absence of superfiuid currents the average equilibrium number M k of Landau quasipar- 
ticles in a state (ka) is given by the element \0'kl^ a kl/ °^ tne ma trix nj^ , 



^-(«> = 5(l-2^F«) (16) 
[compare this expression with the corresponding Eq. @ for normal matter] . 



B. The system of kinetic equations 

To obtain the kinetic equation let us slightly perturb the system. Since our aim in this section is to determine the 
kinetic equation in the linear approximation, we may assume, without any loss of generality, that the perturbation 
varies with coordinate r and time t as e i ( gr ~ u ' t \ where q and w are the perturbation wave vector and frequency, 
respectively. To use the Landau theory of Fermi-liquids we have to assume, in addition, that q <C kFi and uj <C jtXj. 

The only non-zero matrix elements, induced by the perturbation, can be written in a compact form as (see, e.g., 

Sn {t) (n)- A + /> A - fc +-' J / > = m kall 0n kal2 (-i 7 \ 



Ha 21 Ha 22 



Here and below we use the notation 



*-=*-§, *+=* + §• (18) 

The matrix <5n^(q) with the elements <fai^ 5n£\ 22 , defined in Eq. (|17[) . can be interpreted as a small deviation 

(0 

from the equilibrium distribution function n^ , caused by the perturbation. The collisionless kinetic equation for 
6vl£ (q) can be found following the derivation of Ref. 24]. The result is: 

w6n® = <5n« ei l) n - n So® n 5cP - _. (19) 

ka ka k^aO k-aO ka k-aO ka ka fc + erO V * 

In this equation 5e£ is a matrix, describing local deviation of the quasiparticle energy from its equilibrium value 
Cu n« It is the sum of two terms, 

ka U 7 



fc2 = A2 + Ai2, (20) 



where the term describes the change of the quasiparticle energy with the distribution function, 

A (0 = ( E*vj f ij (k,k') 6n^, n £ k , V^(k,k') ^« 12 ^ ( 

and the term is responsible for the interaction of quasiparticles with the self- consistent electromagnetic field, 



Here V and A are the scalar and vector electromagnetic potentials, respectively. It is assumed that in the unperturbed 
system V = and A — 0. 
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Equations (fT9|) -(l22 |) are trivial generalizations, to the case of superfluid mixtures, of the corres pon ding equations 
presented in Ref. [24[. The only non- trivial point is the expression for the coefficient Q!j. In Ref. 24] a, is equal to 
the quasiparticle electric charge ej. This result is valid only for a one-component Fermi-liquid. As is demonstrated in 
the Appendix A, for a mixture of Fermi-liquids one should instead write 



_ mi m-^ 



(23) 



where is the relativistic entrainment matrix at zero temperature, given by [48l . |4£ 



Gi 



(24) 



In this equation &y is the Kronecker symbol and the matrix G{j is defined in Eq. 

The kinetic equation (|19l) consists of four coupled integral equations. Their solution determines the matrix 5nl (?) 



[i.e., the four functions S; 



Sn {€> 1 



The system (|T9|) can be substantially simplified by introducing a new set of variables 



Sri 



(i) 
ka+ 



"feu 11 
0n kall 



Sn 



W 

fccr 22' 
l fccr22' 



r5s 



(0 = 

fccr+ 



(i) 
7, *<t12 



1 W 
'fco-21 



<5s 



(0 

ka- 



5n 



(i) 

ka 12 



8n 



(i) 

ka 21 



(25) 
(26) 
(27) 
(28) 



with the obvious symmetry properties [see Eq. (|I7I) ] 



5?( 



fc(T+ 



5s 



l ka- 

W 

fc(T+ 

(*) 
S fc(T- 



-5n 



-fc-IT +' 



= ft 



l -fc-c 



5s 



-fc-CT+' 



= 5, 



'-fc-c 



(29) 
(30) 
(31) 
(32) 



Using these variables, the system of kinetic equations (|I9p can be rewritten in the form 



x W 




/C-j_ ^K-j- 


p(«)t(<) 
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+ Ai 
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- ^fc„ 




U 5n ka- 






MO 


p(«)t(«) 
^fc + ?fc+ 


^fccr 



(i) 
k + 



2 A, 5s 



(0 

fccr- 



A,; 
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F, 



(i) 



(i) 
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?fc_ 



?fc_ 



(0 

fccr- 



w (5s 



(i) 

fccr — 



A; 



F, 



fc_ 



(i) 



(i) 



fa 



(i) 
ka+ 



f 2 Ai 5n 



(i) 



ft? (i) 
on fccr+ 



- r fc+'Sfc + 



r k- r k + 



A 



(i) 
fccr 



?fc_ ?fc + 



5s 



fccr — 



^fc_"5fc_ "f" ^fc + < 5fc + 



(i) 

fccr+i 



o. 



(0 

fc(T-' 



(33) 



(34) 



(35) 



(36) 



where O fc j + and O fe 



equal 



fccr+ 



E 



V^(k,k')5s^ + , 



o 



ka — 



(37) 
(38) 
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and the functions and are defined as 

Ka Ka 

V£ = 2e 4 y+ £ r(k,k')Sn%l_, (39) 

k'a'j 

41 = 2a^-^/-(fc,fc')^ + . (40) 

The system (l33)) - (j36|) contains all information to calculate linear gauge-invariant response of the two-component 
superfluid Fermi-liquid to a vector (e.g., electromagnetic) perturbation. 



C. The particle current density 



Let us sum Eq. over k and a. Then, using Eq. (fT3")) one obtains that, with the accuracy up to quadratic terms 
in q/kpi <C 1, the last two terms in the right-hand side of Eq. (|3~4")l cancel out and we are left with 



ka ka 



fc + -' v fc_0 



A 



(0 

ka 



E[€-$ 



fccr 



Here we also made use of Eq. (|16l) . To proceed further, we expand £^ in Taylor series 



8n 



(41) 



g 

fc F * 



(42) 



where is the velocity of Landau quasiparticles on the Fermi surface, Vi = VFi (k/k). 

Eq. (|41|) can be additionally simplified with the help of Eq. (|40|) and the expansion (1421) . Neglecting all terms of 
the second and higher orders in q/k-pi, one obtains particle conservation law 



LuSrii = qj i 



with 



Sni = \ E 



ka 



Ji = 



Je 



Ok 



k'cr'j 



rii 

- an — A. 

rrii 



(43) 

(44) 
(45) 



Because N^o lB * ne isotropic equilibrium distribution function, one has dAf^/dk — [dN^l/dk] (k/k), so that Eq. 
can be rewritten as 



Ji = 



2 ^ 

ka 



s$U - E ^ # ^ ^ 



Ok 



fc/ k'a' 



a t — A. 

rrii 



(46) 



where the second term in square brackets was integrated by parts and we made use of the expansion © . 

The quantities 5rii and j i in Eqs. (|43|) -(j46 j) can be interpreted as the number density perturbation and particle 
current density, respectively. To prove that this is really so, we employ Eqs. (f2"5j) and (|26|) and the symmetry relations 
(|29]l and ([30|). and rewrite Eqs. (J44J) and (J45J) in the familiar form, 

(47) 

a^A. (48) 



6n t = SnQ u , 



ka 



E 

ka 



feVj 
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One sees that 5ni is indeed the number density perturbation. Thus, as follows from the continuity equation (|43l) . j i 
is the particle current density. Notice that, j i is formally given by the same expression as for normal (nonsuperfluid) 
Fermi-liquid (see, e.g., Refs. [H, Ef|). For a one-component Fermi-liquid this was first demonstrated by Leggett [H ]. 

Since the quantities Srii and j i are observables, they must be invariant under gauge transformations. This property 
of Sni and j i is very important for the consideration below. For a one-component Fermi-liquid the gauge invariance 
of Sni and j iy defined in Eqs. (|44l) -(j48 |) . was explicitly demonstrated by Betbeder-Matibet and Nozieres [24[ [see their 
Eq. (28)]. 



III. LANDAU FERMI-LIQUID EFFECTS AND GAUGE INVARIANCE 



In this section we study the response of the system to an applied harmonic electromagnetic field. Our aim will 
be to express the exact polarization tensor for a strongly interacting two-component Fermi-liquid through that of 
noninteracting liquid [for which the Landau quasiparticle interaction / y (fc, k ) = 0]. 



A. Gauge-invariant expressions for Sni and j i 



Assume that V and A are some self- consistent (not external) vector and scalar electromagnetic potentials in an 
arbitrary gauge. The vector potential A can be presented in the form 

A = A\+ A tT , (49) 

where A\ = q (qA) / q 2 is the longitudinal component of A directed along q and A tI = A — A\ is the transverse 
component. 

Performing a gauge transformation 

V = V + ^ =V -iuxj), (50) 
at 

A = A-V(f) = A~iq(f> 7 (51) 

one can choose the new potentials V and A in such a way that A\ = (i.e. (f> — —iA\/q). In this new gauge the 
relation between Sni, ji and V , A can generally be written as 

8m = P$ V, (52) 

U = PW + l p(i) (53) 

where Pj and P t ^ are, respectively, the exact longitudinal and transverse polarization functions for a strongly 
interacting Fermi-mixture. The function pj^j is related to Pj by the continuity equation ((43]), 

if = ^ (54) 

A direct (but not easy) way to obtain the quantities Pqq , P/*\ and P t ^ is to solve the system of integral kinetic 
equations (1551) - (f56]l for Sn^_ and Sn^,, and then to make use of Eqs. (pH)) and (1431) . 
Employing Eqs. (|50|) and (|5T|) . one can rewrite Eqs. (|52|) and (|53]l in the original gauge, 

<5n, = pW ^ - ^ , (55) 

U - A tI + I P« - | A,) . (56) 

Clearly, these expressions are gauge-invariant. The knowledge of P 00 , P, v , and P tr allows one to determine the 
longitudinal ei and transverse £t r dielectric functions of the liquid, 

i i 
£ tr = 1-^ E 6 ^- (58) 
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For noninteracting Fermi-liquid, for which k ) = 0, Eqs. (l55t and (|56|) take the form 

Sn t = 11$ f V - ^ , (59) 

U = n£U tr + £n»(V-^A (60) 

where the exact polarization functions of noninteracting matter are denoted as ITqq , nj J , and Il| r . These quantities 
were carefully analyzed in the literature (see, e.g., Refs. 0, [50l - l52j ); they depend on a number of parameters, in 
particular, on the particle Fermi- momentum p-pi and on the mass nil. 



B. Landau quasiparticle interaction and the polarization functions 

(i) (i) (i) 

Let us calculate the exact polarization functions P 00 , , and P tr under the simplified assumption that all the 
Landau parameters except for and f\ J vanish, ff J — for I > 2. Eqs. and (|4"0")) can then be rewritten as 

y W = 2e % V + 2 £ $ 5n 3 = 2e t V$ , (61) 

3 

41 = 2a, — - E lAr # E *' - 2 ei ^S, (62) 

^— ^ /cfi kfi , to, 

J J fc'cr' 1 

where we employed Eqs. (f29|) . (f30|). (|44|. and introduced the effective scalar IQj and vector A^g electromagnetic 
potentials. One can express them in terms of real electromagnetic potentials V and A and gauge-invariant observables 
Sni and j i by making use of Eqs. (j44|) and (|45|) . respectively. The result is: 



y(i) _ 



^+ 1 E/o J ^' ( 63 ) 

J 

A< & = A ■ tE^jV (64) 

To obtain Eq. (|64l) we use Eq. (|7|) for the effective mass to* and the expression (|23[) for ct^ . The detailed derivation of 
Eq. (64) for a one-component Fermi- liquid is given in Appendix B. The matrix 7^ in Eq. (|64[) depends on the Landau 
parameters f\ 3 and equals 

o 

m 

lit = "tt" (GuGij m + m,j - Gu nj) , (65) 
lij = „ J G tj (G u m l + dj rrij - n t ) , (66) 

with 

Si = (Gu rrii + Gij rUj - m) (Gij mi n { + G tj m 3 rij - m rij) . (67) 

Here the symmetric matrix Gy is defined by Eq. (J9)). In Eqs. (|65l) - (f67|) the indices i and j belong to different particle 
species, i ^ j. For instance, if i = 1 then j ' = 2 and vice versa. 

Now, let us return to Eqs. (|61l) and (1()2"1) . The quantities VS and A^ a , entering the kinetic equations ([33]) -([55]), 
depend on the effective potentials and A^g in exactly the same way as in the absence of Landau quasiparticle 
interaction [when f l i(k,k') — 0]. The only difference is that in all equations one should replace the bare mass rrii 
with the effective mass to*, or, equivalently, replace nij with p-pi/v-pi- Moreover, as follows from Eqs. (|6H)) and (|64l) . 

the gauge transformation properties of V^ s and A^g coincide with that of, respectively, V and A [see Eqs. ([50| and 

(|51|)]. Consequently, the relation between Sni, j i and V^g , A^g in an arbitrary gauge is given by the same expressions 
as for a noninteracting Fermi- liquid [compare with Eqs. (|59l) and (1601) ]. 

Sr* = n« V (68) 



<1 
Q_ 

lo ' \ q 



nw^tr + ^nwfvw-^). (69) 
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Let us emphasize once again that here the quantities IIqq , , and 11^ should be understood as the functions of ppt 
and m* (or ppi and vpi = ppi/m*) rather than the functions of p-pi and m^. These two sets of variables are equivalent 
only for noninteracting Fermi-liquid, when m* — mi. 

Now we are ready to find the polarization functions Pqo , P\ ■> an d P^ f° r a two-component, strongly interacting 
Fermi-liquid. For that, we compare the general expressions (|55[) and f|56[) for Srii and j { with the corresponding 
equations (|6"5)) and Taking into account, that and A^2 are given by Eqs. (155)) and we obtain 



/>(') = 1 . . . —L (70) 

00 p-p- -p- v-n (3) -p- Y- n (l) - y ■ v-n (l) n (i) + v - v n (l) n (i) ' 

cjPj Pj Xjj ^OO e J A" 1J -00 AiJ Aji 11 00 11 00 ' W AjJ 1J -00 11 00 



n Q ^Pj + ej Xij U-qq Pi n, 



pM - ^ : : : l n\\ 



while P, can be found from Eq. (I5H) . Here Xij = /o' 7 — {u 2 /q 2 ) Jij and the indices i and J refer to different particle 
species, I ^ j. Eqs. (|70|) and (|7T|) will be analyzed in a subsequent section. 

Here it is convenient to make a few comments, concerning the scheme of calculation of the polarization functions, 
suggested above. First, in principle it would be possible to extend this scheme taking into account the Landau 
parameters j\ 3 with I > 2. However, for that one needs to know how the 'noninteracting' system [with f' lJ (k,k) = 0] 
responds to a general perturbation (not just to the electromagnetic field). This is a complex problem without any 
simple solution, such as in the case of Z < 2. 

Second, the approach described above can be easily generalized to calculate the axial-vector response of the system. 
For that one needs to introduce the spin-dependent part of the Landau quasiparticle interaction in the kinetic equation 
(I19|) . For a one-component superfluid Fermi- liquid, the axial- vector polarization functions were calculated by Leinson 
Q under the simplified assumption that only nonzero Landau parameters are go and g\ (see Ref. [t| for the definition 
of go and gi). The generalization of his results to the case of mixtures is straightforward. The more difficult problem 
would be to estimate the effect of tensor quasinucleon interactions on the axial- vector response of the system. We 
hope to address this problem in our subsequent publication. 



C. Various limiting cases for Pqq and P$ 

We consider first a one- component Fermi-liquid. The particle species indices can then be suppressed. The polariza- 
tion function P o for a one-component Fermi- liquid was studied by Leggett [23| in the limit of qv? <C A and ui <C A. 
Both polarization functions Poo and Pt r were studied in the recent paper by Leinson Q at arbitrary qvF <C and 
lj <C /i. In his analysis, Leinson took into account only the harmonic Vo of the pairing interaction and assumed that 
Vi = for / > 0. Our results will be compared with the results of these two authors. 

For a one-component Fermi-liquid Eqs. (|70|) and (fTTj) are essentially simplified 

^oo = (72) 



p - x IIoo 

en tr 

e — 7 n tr ' 



P< = (73) 



with 



171771 £ in A S 

7 = -ns-/i» ( 74 ) 

ftp 

X = To r 7- (75) 

q 

To obtain Eqs. (|T2"j) and ([75)1 we put = and f\ 3 = for i ^ j in Eqs. (1701) and (J7TJ) and then suppress the particle 
species indices. 

Eq. ([72")) coincides with Eq. (68) of Leggett [23[ and with Eq. (63) of Leinson Q in the limit of qvF <C A and 
u <C A. In the other limiting case, when qvF <C A and u> > 2 A, we reproduce the result of Leinson Q, see his 
Eq. (82). However, at arbitrary qvF and ui our Eq. (fT2")l differs from the general equation (55) for Poo, suggested by 
Leinson Q. 
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His Eq. (55) depends on two complicated combinations of integrals, Q and P, that can be easily expressed through 
each other only in the limits discussed above. In contrast, our Eq. (|72p depends solely on the polarization function 
Hoo(pf,Vf) of noninteracting one-component Fermi-liquid, which is the direct consequence of gauge invariance of the 
quantities 5n and j. 

Now let us examine Eq. (j73[) for P tr . Knowledge of P tr allows one to calculate the transverse-current autocorrelation 
function Kt, 

X T = — + -• (76) 
e m 

This quantity was derived by Leinson (see Eq. (86) of Ref. j9j). Notice that his Eq. (86) contains a misprint (53|; one 
should replace Vp with pp/m in this equation. After correcting the misprint, Eq. (86) of Leinson coincides with our 
Eq. Gii. 

To further check Eq. (|73|) . we look at the static limit of Pt r , assuming that u = and qvp «C A. In the static 
limit the particle current density j is generated solely by a motion of superfluid liquid component (i.e., the normal 
component is at rest). It is given by [43] 

j = — V s , (77) 
m 

where p s is the superfluid density and V s is the superfluid velocity. The velocity V s depends on the gauge-invariant 
combination of the phase (p of the Cooper-pair condensate wave function and on the electromagnetic potential A (see, 
e.g., Ref. [47j and Sec. IV for more details), 

V s = ^-{V^-2eA). (78) 
Ira 

In the transverse gauge in which A\ = (qA) /q = 0, the phase ip = 0, because it can only depend on the scalar (qA) = 0. 
Then it follows from Eqs. ([77]) and (J75J) that 

j = A tr . (79) 

For a one-component noninteracting Fermi-liquid 

ft = mn[l-$(T)], (80) 

where $(T) is a function of temperature [for more details see, e.g., Ref. [44j |. where this function was denoted by 
f(T)}. Using Eqs. flM}, and (JHOJ), one obtains 

„ / * e n vp (1 — $) , , 

n tr (pf,vf) = £-i ^. 8i 

This equation together with Eqs. (O and ([73]) gives 

Ptr - "m(l + Fx$/3) ! (82) 

where Pi = (m*pp/'K 2 ) f\. Comparing Eqs. (l79l) and (|82|). one can determine the superfluid density p s for interacting 
Fermi-liquid, 

m 2 «[1-$(T)] 

Ps " T tr " i + F!$(r)/3- (83) 

It coincides with the corresponding Eq. (72) of Leggett [Hj]. (Notice that, Leggett calculated the 'normal' density 
p n = mn — p s .) Thus, we demonstrate that in the static limit our Eq. (|73l) reproduces the well-known result of Ref. 
44]. 

We now turn our attention to the two-component Fermi-liquid and discuss first the static limit (u) — and qvpi 
Aj). To obtain in this limit one should use Eq. dHJ with U^J given by expression, similar to Eq. (|8ip. 

n r (PFi,WFi) = , (84) 

PFi 
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where $i is the same function of temperature as $ (for more details, see Ref. [38j])- 

On the other hand, in the static limit the general hydrodynamics of superfiuid mixtures is applicable, which states 
that in the absence of normal current (see, e.g., Ref. |36j|). 



^ rrn 



(85) 



Here pij is the entrainment matrix (also termed the Andreev-Bashkin or mass-density matrix) and V a i is the superfiuid 
velocity for i-th particle species. It can be expressed through the phase ipi by an equation, similar to Eq. (|78p . 



V si = -!- {V<fi - 2eiA). 
2m; 



(86) 



Again, as for a one-component liquid, in the transverse gauge ipi = and from Eqs. (|55|) and (|56")l we have 



J i 



&i Pii . &j Pij 
2 1 I 

mf m, m,4 



(87) 



or, in view of Eq. (|56[) . 



P (0 _ 

^tr — 



e» p v 



+ 



e j Pij 

mi rrij 



(88) 



In Eqs. (|87|) and (|88|) indices i and j refer to different particle species, i ^ j. Comparing Eq. (|88l) with Eq. (|7T|) and 
taking into account Eq. (|84p , one can determine the entrainment matrix p,j and verify that it coincides with the result 
of Ref. (HI, obtained in a quite different way. (It should be noted that in Ref. [38] it is additionally demonstrated 
that the higher harmonics /, with I > 2 do not contribute to pij.) 

The next interesting limiting case is realized if one particle species, say, i = 1, is charged while the other is not 
(e 2 = 0). It follows then from general equation ((701 that 



P, 



(i) _ 



III! 



eill$ 



1 



X22ll$ 



d(2) _ 
00 — 



X22 -li _ Xll ^00 ~ Xl2 X21 1%) n 00 

„ v n( 1 )fr(2) 
ei X21 Hqo n oo 



Xll X22 11 11 00 



X22 11 



Xn Hqo 



Xl2 X21 tl li 00 



Xll X22 11 11 00 



(89) 
(90) 



Here the function n Q is independent of the electric charge ej, n g = n^g/e^. An analogous expressions for P t 
and P tl . can be obtained with the help of Eq. (|71[) . One sees from Eqs. (l89l) and (|90| [and from the corresponding 



equations for P^- 1 ] that neutral particles not only modify the polarization function PgQ ; tr of charged particles, but 

(2) 

also respond themselves to electromagnetic field, because P 00 tr ^ 0. However, the neutral particles do not contribute 
to the dielectric functions of the liquid, since, as follows from Eqs. (|57|) and (|58|) . e\ and e tr are given by 



(i) 



Etr 



,(1) 



4-7T 

1- — ei p t w. 



(91) 
(92) 



Finally, we mention another interesting property of the general solution ([70)1 and (fTTj) . Assume that the polarization 



functions H 

of n w n (j) 

oi ii 00 ii 00 



(i=l,2) 
30 

[or H 



[or n; 

tr lx tr . 



(i=l,2)i 



of noninteracting Fermi-liquid are small so that one can neglect the terms of the order 



. Then, as follows from Eqs. (1701) and (fTTj) . in the linear approximation, 



P 



« 



oo 



H 



« 
oo 



or P 



n 



(i) 



(93) 



That is, the polarization functions are not modified by the first two harmonics and f[ J of the Landau quasiparticle 
interaction. In the next section we demonstrate that in some cases this conclusion is correct even if we take into 
account all harmonics f, 3 with I > 0. 
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D. Transverse polarization function in the Pippard limit 



In the Pippard limit we have qvpj 3> Aj, qvpj 3> w, and qv-pj -C fJ-i- It is especially important to know the 



polarization functions Pqq and P^' in this limit because they are required, for instance, for calculating the kinetic 
coefficients of a multi-fluid Fermi-mixture in neutron-star cores (33l . l34j . 

In the first approximation, Pqq does not depend on the gap Aj and the frequency w and agrees with the corre- 
sponding expression for normal matter, describing the ordinary static screening of particles. It can be easily obtained 
from Eq. (|T0[) if we notice that for a normal one-component noninteracting Fermi- liquid one has (see, e.g., Ref. [45|) 



n 



(<) _ 

00 — 



e iPFi 
7T 2 W Fl 



(94) 



Strictly speaking, Eq. (|70p that we employ, was derived under the assumption that the only harmonics f and of 
the Landau quasiparticle interaction are nonzero. However, one can easily verify that, for a normal Fermi-liquid, P 00 
is still given (at small q and lu) by Eq. ([701) even if we allow for higher harmonics j\ 3 with I > 2. For a one-component 
Fermi-liquid this was demonstrated, for example, in Ref. pF 



Now let us consider the transverse polarization function P^/ . In the Pippard limit the polarization function is 

small, IL^ = 0[Ai / (qvpi) + uj / (qv-pi)]. It follows then, from the discussion at the end of the preceding section, that 

the first two harmonics /q 7 and of the Landau quasiparticle interaction have no influence on P t */ , so that it is 
given by Eq. (|9"5|). 

Below we demonstrate that this result remains correct even if we take into account other harmonics fl with 
I > 2. In the static limit (cj = 0) this was first shown by Leggett (lij . More precisely, we prove that the transverse 
polarization function for a noninteracting system coincides with the function for a system with an arbitrary harmonic 
fi 3 {I > 2) switched on. 

We consider first the simplified situation in which the pairing potential V^'(k,k ) is a constant; the generalization 
of our results to the case of an arbitrary (k, k') is briefly discussed at the end of the present section. In other 
words, we take into account only the first term Vg^ in the expansion (1121) of V^\k,k') in Legendre polynomials and 
neglect all other terms, V t =0 for I > 1. In that case the functions O^' and O ka _ in Eqs. (|33|) - (|3ol) are some 
constants depending on the scalars qA and V (see, e.g., Ref. 47]). In the transverse gauge, where qA = qA tI = 0, 

= 0. The solution to the system 
(q,u, Aj) 



r(i) 



and in the absence of scalar electromagnetic potential V, they vanish, Oj£, = Oj^._ 



of equations (|33|) - (l36l) is then simplified and we have for Sn^l , (q, to, Aj) and 5n£ 



(0 



Ai 



5n ( *l_(q,uj, A, 



■«v, 



(i) 



Mi 



(i) 



Mi 
IT 

Ma 
D ' 



(95) 
(96) 



where 



Mi 



1 



Mo = - 



2 

Ai) 

1 
2 



4A 2 



F 



(i) 



-4(ffu,) a ^ i)2 + (fl»i) a w 



F 



2 , .2 



9»< C k 



F, 



(i) 



D = IE: 



E. 



0) 



4^ 2 -u 
F^ - a; - £0 ) ( F, 



F, 



« 



F 



iO F 



p(0 
fe_ 

(i) 



F, 



f[ 4) - uj - iO 



F 



(i) 



SO 



(97) 



(98) 



(99) 



In Eqs. dSU) and dHSJ) V$ and Aj£ are the smooth functions of A;, defined in Eqs. (|59")l and P0|) . respectively. For our 
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problem they can be rewritten as 



k'a'j 



feVj 



(100) 
(101) 



where 9 is the angle between k and A:' and I > 2. 

As already mentioned above, Leggett [HJ showed that the static function pfy(q, 0, A,) is not affected by the Landau 
quasiparticle interaction. Therefore, it is sufficient to analyze the difference P^*\q, to, A*) — P t "? (g, 0, Aj) and prove 
that it is independent of /, u . Using Eq. d46l) , we obtain 



(102) 



ji{q, u, A,) - j t (q, 0, Aj) = [p« (g, w, A,) - P« (g, 0, Aj) 

= 5 E^: h^+^^^-^+^o.A,) 



The function fcj^) , (g, w, Aj) — (Jnj^, (g, 0, A,) is nonzero only in a narrow region near the Fermi surface, when k ~ fepj. 
Furthermore, because of the denominator D [see Eq. ([55)1]. this function has a sharp maximum at guj < (w + Aj). 
Introducing the longitudinal k\ || q and transverse fc tr -Lg components of the vector k = k\ + k tr , this inequality can 
be rewritten as 



rrii{ui + Aj) 



(103) 



The main contribution to the integral (|102p comes from k\ satisfying Eq. (I103p . Keeping this in mind, it is straight- 
forward to verify (see, e.g., Ref. [5l|) that for noninteracting Fermi- liquid one has 



j i (q,u,A i )-j i (q,0,A i )=O 



A ? ; 



(104) 



Now let us analyze whether the Landau quasiparticle interaction influences this result. For this purpose we inspect 
the terms in the function Sn^Aq, uj, Aj) — Sn^Aq, 0, Aj) which depend on f^. From Eqs. ([95]). (llOOp . and (jlOlf) it 
follows that they have the form, 



k' a' 



kk' 

kFikFj 



6n%l + and /J = £/«P 



kk' 



8n 



(J) 

k' a' 



(105) 



By demonstrating that these integrals are quadratic in (cj + Aj) / (qvYi) we prove that P$ — Il^f {pFi, v-pi) in the 
Pippard limit. Below we consider in detail the first term in Eq. (|105l) ; the analysis of the second term is similar. 
The first term can be presented as 



k'cr' 



k' a' 



k' a' 



kk' 

kFikFj 

kk' 

kpikpj 

kk' 

kFikFj 



Sn^Jq^A^-Sn^JqA^) 



Sn k i, + (q,0,0). 



(106) 



One can easily verify that the last integral here vanishes at I ^ 1 and thus can be omitted. Furthermore, because 
of the same reasons as those discussed below Eq. (|102[) . the main contribution to the first two integrals in Eq. (|106|) 
comes from such a region of k' , that 



K<^±^«K^k ¥j . 



(107) 
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From symmetry arguments it follows that the functions in square brackets in these integrals can generally be written 
in the form 

[...] = (k'A) G (k'q, k') , (108) 

where G{k'q 1 k' 2 ) is a scalar function. 

Now, using Eqs. (|T03| . (fT07|) . and (fT08)l , one may write 

E ff a (^t) [• ■ J - E /f « G "0 
= E /? p - f rf) g «», *') + o ( ( " + A - )( " + A '' 

1 tr ^ 

= offr + AQfr + AjA (109) 
V qv F i qv Fj J 

Here we used the fact that the integral over directions of k' tT vanishes for I ^ 1. Thus, we demonstrate that / ~ 

(lu + Ai)/(qvFi) x (u> + Aj)/(qvFj) and, consequently, the Landau quasiparticle interaction has no influence on 
in the Pippard limit. 

(i) 

The consideration of this section is simplified since we take into account only the harmonic Vq of the pairing 

potential. However, it seems plausible (though we have not checked it in detail) that the inclusion of other harmonics 

(i) 

Vy with I > 1 will not change the result. In principle, a prove of this more general statement should be similar to 
the prove presented here, but the equations to be analyzed, are much more complicated. In particular, the function 
5t%1 + wu 1 depend, in addition, on the integrals O^, and Oj^_ [see Eqs. (|57|) and (|3"8"])]. 

Finally, let us make a comment concerning the coefficients and of the first harmonic. From the analysis 
presented above it is clear that the coefficient f\ J plays a very special role, because the integrals in Eqs. (|106[) and 
(I109[) do not vanish at I = 1. It might think that the situation with the coefficient of the pairing potential is 
analogous, so that the integrals 

Y, VP Pi ( cos 9)S8$ a Jq,u, A t ) and ^V^ 3 P^cosfl) Sa^ a+ {q, u, A*) (110) 

k'a k'cr 

are nonzero. However, this is not the case; they vanish due to the symmetry relations (|31[) and (1321) . 



IV. A PHENOMENOLOGICAL APPROACH TO KINETIC EQUATION AT SMALL q AND uj 

In this section we analyze the kinetic equation in the 'quasiclassical' limit, qvpi <€. Aj and w Aj (see, e.g., 
Refs. [lillil])- This limit is especially important for various applications, for instance, to study low- frequency long- 
wavelength collective modes propagating in superfluid matter and to calculate kinetic coefficients. 

For a one-com pon ent Fermi-liquid the quasiclassical limit of the matrix kinetic equation (|19[) was thoroughly 
examined in Ref. [24j . It was demonstrated, that Eq. (|19[) can be substantially simplified by introducing a concept of 
Bogoliubov excitations. In particular, the kinetic equation for Bogoliubov excitations acquires a scalar (rather than 
matrix) form. For a two-component Fermi-mixture the analysis is quite similar. Here we do not attempt to perform 
it; an interested reader is referred to Ref. [24] for more details. Instead, we follow a more intuitive phenomenological 
approach allowing to formulate the kinetic equation in the non-linear regime (in contrast to Ref. [24j . where the kinetic 
equation is studied only in the linear approximation). We call this regime non-linear in a sense that, for instance, it 
allows us to study the nonequilibrium variations of the energy gap which are comparable to Aj. However, because we 
use the Landau theory of Fermi-liquids, we still assume that the quasiparticle distribution only slightly differs (in the 
vicinity of the Fermi surface) from the step function. Our results will be discussed and compared with those available 
in the literature in the end of the present section. 

In the previous sections the quasiparticle momentum was denoted as k. In the Appendix A we demonstrate that 
in the presence of electromagnetic field k is actually a generalized momentum. To distinguish between k and the real 
momentum, the latter will be denoted by p. It is more convenient to use p instead of k in the consideration below. 
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A. Local analysis 

In the quasiclassical limit one can assume that a Landau quasiparticle (or a Bogoliubov excitation) with a certain 
momentum p possesses, at the same time, a certain coordinate r. Consequently, such quantities as the distribution 
of Landau quasiparticles (Bogoliubov excitations) , their energy, or the energy gap can be considered as functions of 
p and r. To find how these quantities are related to each other it is sufficient to analyze the system locally. 

Let us consider a two-component Fermi-liquid out of thermodynamic equilibrium. To simplify the problem we neglect 
for a while the electromagnetic field, assuming that the liquid is composed of neutral particles. The electromagnetic 
effects will be taken into account in Sec. IVC. Our aim will be to calculate the energy density E of supcrfluid matter 
in the neighborhood of a point r. In the vicinity of this point the matter is almost homogeneous. Thus, it can be 
approximately described by a uniform Hamiltonian H (see, e.g., Refs. [38l. 144 149|). 

H-J2 fa N i = #LF + ^pairing, (HI) 

i 

where Ni is the number density operator; /tj is the nonequilibrium analogue of the chemical potential /ij to be 
determined below; Hlf is the Fermi-liquid Hamiltonian, 



pert 

+ \ £ f ij M{^Hl-e^)(a^l-9^), (112) 

pp' (Td'ij 

and i/pairing is the pairing Hamiltonian. In the presence of superfluid currents it is given by (see, e.g., Ref. [38]) 

-^pairing — V Q % \PiP ) "p'+Qit -P'+Qii ~P+Q t i P+Q.t' \ 11C >} 

pp'i 

Here, 2Q { — 2miV S i is the momentum of a Cooper-pair in the condensate. It is related by Eqs. (|143l) and (11451) to 
the quantity introduced in Sec. IIB. The matrix element Vq (p,p') in Eq. (II 1 3|) describes scattering of a pair 

of Landau quasiparticles from the states (p + Q 4) t), (— P + Qi, I) to states (p' + Q i: f), (— p' + Q i7 1). In Ref. [H[ it 
is argued that Vq (p,p') « V^ 4 )(p,p') [up to small terms ~ (Qi/fcpj) 2 ]. 

To find the energy density E one needs to diagonalize the Hamiltonian (|111[) . To do that, we rewrite Eq. (Illlj) in 
terms of Bogoliubov operators bpl, defined as (see, e.g., Ref. (38j ) 

p+Qit ~ v °p+Q^ + v p °-p+QiV > 

„W - v {i) d<) _ ,,(») I,Wt hie) 

p+Qii — p p+Q,l p -p+Qit' ^ > 



where u p ^ and i;^ are even functions of p, 



normalized by the condition 



u$=u% v?=v% (116) 



u^ 2 +v p ^ = l. (117) 

(i) (i) 

Generally, the coefficients u p and v p are complex. However, at some moment of time they can be chosen to be real 
in the vicinity of our point r by a suitable phase transformation. 

Being expressed through the Bogoliubov operators, the Hamiltonian (|111[) takes the diagonal form. Thus, one gets 
the following expression for the energy density 



e -£m^=£ [4+QiO 

i peri 

pp'aa'ij 



i) _ f)(i) 

p+Q, °p+Q, 



r (p + Q i)P ' + Qj) (« p % Qi - K'U - 



+J2 v(l) M 4HH' v p' 

' ->P Q -'"p q ) >P Q ^ +0 J 
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Here Fp\.Q. is the distribution function for Bogoliubov excitations with momentum (p + Qi), 



F, 



- /ii,Wt i,W 



p+Q, — M "p+Oit"p+Oi1"^ = 
while A/p^g is the average number of Landau quasiparticles in a state (p + Q i ,<j), 

r(0 - /i„(0t „« i\ _ /i„Wt „W 



(119) 



Jsl p+Q, ~ " a P+Qit a P+QiV' ~ " a P+Qil a P+QiV' 

= v^+^T^ Qi -v^T% +Q . (120) 

If we were in thermodynamic equilibrium we could easily find the unknown functions Fp\.Q and Up in Eq. (|118p 



by requiring minimum of the free energy F, F 



p+Q . ) "p 



E 1 — jUini — /ii2^2 — rS*, where the entropy 5[^ p+ q.] 
is the functional of only F^q (see Ref. [38| for more details) . Since we are not in thermodynamic equilibrium, the 

(i) 

distribution function for Bogoliubov excitations F p ^Q in our local analysis should be considered as a given 'input 
parameter' [it can be found from the corresponding Boltzmann kinetic equation (|132[) , see Sec. IVB]. To determine 
Up \ we assume that even out of equilibrium F still has a minimum as the functional of Up (at fixed ). This 

assumption, though plausible, cannot be proven in our phenomenological approach. However, its validity can be 
justified by comparison with the results of the strict microscopic theory (see Sec. IVC). 
One obtains from the minimization procedure 



1 



U ^=2 11 



H, 



p+Q, 



p+Q, 



2<£ 



p+Q, 



H 



(i) 

p+Q, 



-H. 



(i) 



p+Q, 



where 



H {i) = £ (i) KfU) 
P P P 



see Eq. ((T|) for the definition of Sp , and 



*(*) 
~p+Qi 



p+Q, 



-H 



-p+Q, 



H. 



p+Qi 



H 



(i) 

-p+Q, 



(i)2 



(121) 



(122) 



(123) 



is the energy of a Bogoliubov excitation with momentum (p + QJ. To verify that &p\.Q is indeed the energy, it is 
sufficient to notice that it is given by the variational derivative of the functional (E — pLiTii) with respect to F p ^Q , 

S(E- 



-p+Q; 



SF. 



p+Q; 



Finally, 2?p in Eq. f| 123[) is the nonequilibrium energy gap. It is defined by the equation, 



2>W = - £ (P,p') (l - - J" 

p' 



(<) 
p'+Q, 



(124) 



(125) 



It can be demonstrated that, in the linear approximation of Sec. II, the quantity ST>p^ = T>p — Aj is related to the 



integral O p i] a+ [see Eq. (J37J)], 



W 



(126) 



Using Eq. (|120[) one can determine the nonequilibrium chemical potential pn from the requirement that the number 
density rii is given by the sum over all occupied quasiparticle states, 



p+Qi 



(127) 



pa 
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The quantities ■A/'I+q.! &p+Q.i an d ^p can be easily found from, respectively, Eqs. (I120p . (|123l) . and (|125[) . once 
le distribution is specified. As sho 

by the standard Fermi-Dirac distribution, 



the distribution Fp\_Q is specified. As shown in Ref. [38j , in thermodynamic equilibrium the function J~p]_Q is given 



^ifO,o = i /T - ( 128 ) 

Here and below the equilibrium function J- p ^ is denoted as J-„\, where tilde indicates that we allow for superfluid 
currents in the system. The function J~ p0 should not be confused with the distribution J~ p q [see Sec. IIA and, in 
particular, Eq. ([5]) for the definition of J~po\- These functions are equal only in the absence of superfluid currents 

Wi = o). 

^io = 4o- (129) 
In this case one also has for the system in thermodynamic equilibrium 

<£« = (130) 
V® = A*. (131) 

The equalities fT5SF |> — ([T5T |> follow from Eqs. Q, ©, and (HDJ) of Sec. IIA and Eqs. l[I2D jl -([I25 jl of the present section 
(see Ref. [38] for a detailed derivation). 



B. Introducing dynamics 

The relations between various nonequilibrium quantities discussed above should be supplemented by the kinetic 
equation for Fp+Q , the continuity equation, and by the 'superfluid' equation, describing the evolution of Q i with 
time t. The kinetic equation for the Bogoliubov excitations takes the standard form, 



dt dp dr dr dp 



The collision integral on the right-hand side of this equation can be easily obtained once the interaction between the 
excitations is known (see, e.g., Refs. [25l - l30l . l54j ) . 



The continuity equation is written as 



^ + divj, = 0, (133) 



where the number density equals 



per pa- 

and the particle current density is given by (see, e.g., Ref. [44 

dH® 



E^^P% 4 - d35) 



pa 



Since A/^q differs from the equilibrium distribution A^q only in the narrow region near the Fermi surface, Eq. (|135[) 
can be linearized and rewritten in the form, similar to Eq. (|48|). or, after some algebra, to Eq. (J46j) . Then, using Eqs. 
(fTTB]) and (fT20l and noticing that 

7 p+Q, -p+Q; ~ ^p+Q, -p+Q s ' ^ ^ 
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Eq. (|135[) can be finally presented as 



h = E Y * 



«i + h E p 



pa 



(137) 



where the matrix Yij is given by Eq. (|24[) . For a one-component Fermi- liquid Eq. (|137l) reduces to the well-known 
expression (see, e.g., Rcfs. 28, 54]), 



n Q P t- 
r 7_v — -^P+Q- 

TO * — ' TO 
per 



(138) 



To obtain this formula we employed Eq. (J7J and the definitions ([9]) and ([241) . 

It is important to emphasize that generally the mass current density rrii j i of «-th particle species is not equal to 
the momentum density Pi, 



(139) 



per per 

However, using Eq. one can check that, due to Galilean invariance of the system, the following equality holds 

E^ii=E p - ( 14 °) 



Now let us discuss the 'superfluid' equation. It has a natural form (see, e.g., Refs. [2J, [28|, [23, [54| , where similar 
equations are written for a one-component liquid) , 



dt 



(141) 



and coincides with the 'superfluid' equation (4.9) of Ref. [28|. (The authors of Ref. [28[ used different notations. In 
particular, our quantity p, is related to their invariant potential $ by fi = [i — <&, where fi is a constant from which 
the authors count energy.) Eq. (| 141[) is also equivalent to the corresponding equation of the Khalatnikov's superfluid 
hydrodynamics [55| . 



dV s / F s 2 \ 



(142) 



To prove this, we notice that the Khalatnikov's chemical potential fj,Kh (per particle mass to) is defined in a reference 
frame in which V s = Q/m = 0, while our potential ji is defined in the laboratory frame. They are related by an 
obvious formula, p, = m (^Kh + V s /2). Thus, Eqs. (|141|) and (|142[) are indeed equivalent. 

Notice that Eq. (|141[) will be automatically satisfied if we express Q { and fa through the wave function phase ipi of 
the Cooper-pair condensate (see, e.g., Refs. [H, H3, H3] ) , 



Qi 

-fa 



1 

2 94 



(143) 
(144) 



One can verify that, in the linear theory considered in Sec. II, the phase ifi is related to the 'zero harmonic' of the 
function Oj£l_ [see Eq. (l38l) ] by the equation 



- — V V (l) Ss {t) 

~ 2A 2.^1 v o 0S ka- 



(145) 



All other harmonics of Ojj. _ are small and can be neglected in the quasiclassical limit [24|. 



19 



C. Inclusion of electromagnetic field and comparison with the previous works 

Eqs. (|TT6]) . (fTTTl) . (H20j)-(IT23l), (|T25]) . (IT271) . (fT28l) . (fT32l) . (IT331) . (IT37I) . (fT39l) . (IT431) . and ([144]) of Sees. IVA and IVB 
fully describe the two-component superfluid neutral Fermi-liquid in the limit of small q and w. The generalization 
of these equations to the case of charged mixtures is straightforward (see, e.g., Refs. (28|, [Hj]). Namely, they remain 
essentially the same if we redefine the quantities Q i and fa [see Eqs. (|143|) and (|144[) ] through the gauge-invariant 
combinations, 

Q l = ^Vpi-etA, (146) 

-^ = lw + eiV - (147) 

Notice, however, that the 'superfluid' equation in the form (|141j) is no longer valid. As follows from Eqs. (|146j) and 
(|147[) . it should be replaced by 

^ = -Vfa + ei E, (148) 
at 

where E = —dA/dt — V V is the self-consistent electric field, which can be found with the help of the Maxwell 
equations. The right-hand sides of Eqs. (jl46[) and (jl4T[) are indeed gauge-invariant since the phase ifi transforms as 
ipi — » ifi — 2ei(f> under the gauge transformation (|50]) and (|5Tj) (see, e.g., Refs. 0, S3|)- 

The equations formulated above in Sec. IV reproduce various limiting cases that were studied in the literature. 
First of all, for a one- component Fermi-liquid our equations coincide with those obtained by Betbeder-Matibet and 
Nozieres [24| and by Aronov and Gurevich [54[ (see also Refs. [28|, HgJ]). Betbeder-Matibet and Nozieres worked in 
the linear approximation but assumed the most general form of the Landau interaction f(k,k') and pairing potential 
V(k,k'). On the contrary, Aronov and Gurevich derived, from the first principles, the fully nonlinear system of 
equations, describing the superfluid Fermi-liquid in the quasiclassical regime. However, they completely neglected the 
Landau interaction [f(k,k') = 0] and took into account only the zero harmonic Vq of the pairing interaction (Vi — 
for I > 0). 

Our equations for a two-component Fermi-liquid were compared only with the results of Sec. II since we did not find 
a discussion in the literature on transport properties of strongly interacting superfluid Fermi-mixtures. We checked 
that in the linear approximation our equations reproduce the quasiclassical limit of kinetic equations (|33p - (|36[) . 



V. SUMMARY 



This paper discusses transport properties of a mixture of two superfluid strongly interacting Fermi-liquids. A typical 
example of such mixture is the matter in the internal layers of neutron stars. To describe the mixture we use the 
Landau theory of Fermi-liquids generalized by Larkin and Migdal [U |43[ and by Leggett [H, [44| to take into account 
the effects of superfluidity. 

Our results are summarized below. 

(i) Working in the linear approximation, we formulate the system (|33[) - (I36|) of kinetic equations, describing the 
collisionless superfluid Fermi mixture in the self-consistent electromagnetic field. To derive these equations, we follow 
the approach of Betbeder-Matibet and Nozieres [24| , who obtained the kinetic equation for a one-component superfluid 
Fermi-liquid. Generally, the system (|33l) - (f3"6")l is a straightforward generalization of the corresponding equations of 



Ref. |24| . However, there is one nontrivial difference concerning the form of interaction of Landau quasiparticles with 



the electromagnetic vector potential [see Eq. (|2"2"|) ]. For a one-component Fermi-liquid a>i in Eq. (|2"2")l is always equal 
to electric charge, ctfj = ej, which is the consequence of the Galilean invariance of the system. On the contrary, for a 
multi-component mixture Ui is generally given by Eq. (|23p . while the Galilean invariance requires only that Eq. ([7]) 
must be satisfied. 

(it) Using the above kinetic equation, we determine the particle current density j i of i-th particle species [see Eqs. 
(|4"5j). (|4l)|) . or (|4"8")) ] and show that it is given by the same expression as for a nonsuperfluid matter. For a one-component 
superfluid Fermi-liquid this was first shown by Leggett [44j . 

(Hi) Assuming that the only two harmonics Jq and of the Landau quasiparticle interaction k') are nonzero, 

we calculate the polarization functions Pqq , P, , and P$ [see Eqs. (|54")l . (1701) . and (1711) 1. and compare them, in various 
limiting cases, with the results available in the literature (see, e.g., Refs. [9j, [23|, |3g, |44j). We demonstrate, that the 
functions Pqo j tr can ^ e expressed through the Landau parameters and polarization functions IIqq 1 tr (ppi,«Fi) of 
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noninteracting Fermi-liquid, for which f^(k,k') — 0. This result is valid for any smooth pairing potential V^(k,k') 
and for all wave vectors q and frequencies u) such that qvpi -C Hj and u) <C fij . 

(iv) We show that the transverse polarization function P$ does not depend on the Landau quasiparticle interaction 

f l j(k,k') in the Pippard limit, when qvpj 3> A,; and qvpj 3> u). In this limit it is given by = Tl[]} (ppi, Dpi). For a 
one-component Fermi- liquid and lu = the same result was obtained previously by Leggett [44] . 

(v) Finally, we formulate a system of nonlinear equations describing the nonequilibrium superfluid mixture in the 
quasiclassical limit (qvpi <C Aj and w< A 3 ). It consists of Eq. (|125[) for a nonequilibrium energy gap, scalar kinetic 
equation (| 1 32|) for a Bogoliubov excitations, continuity equation (|133[) . and of Eq. (|148[) for the superfluid velocity. In 
the linear approximation this system is completely equivalent to kinetic equations formulated in Sec. II and in Ref. 
[2~il ]. Moreover, we verified that it reproduces the nonlinear equations derived from the first principles by Aronov and 
Gurevich 54 1 (see also Refs. [HI, To simplify the problem, these authors neglected the Landau quasiparticle 
interaction and assumed that the pairing potential is a constant. 

The results obtained in this paper can be useful in a variety of applications. For example, the polarization functions 
can be used to study low- frequency (w -CuJ long- wavelength (qvpi <C fJ-j) collective modes in superfluid matter of 
neutron stars [see, e.g, recent papers pll l22j for an example of such studies in the normal matter]. 

Also, the complex parts of the polarization functions Pqq and P t , determine energy losses in the Cooper-pairing 
neutrino emission process. This process regulates thermal evolution of neutron stars 0, [56 58] and is especially 
important in the crust of accreting neutron stars, exhibiting X-ray superbursts [Bl^. l6(il|. 

Next, the kinetic equation derived in Sec. II can be used, with minor modification, to study the axial-vector response 
of a superfluid Fermi-liquid (see Sec. IIIB for more details). This problem is also important in application to the 
Cooper-pairing neutrino emission process 0, H, H, [t| . 

One needs the polarization functions in the Pippard limit since they describe the plasma screening of the interaction 
between charged particles (e.g., protons and electrons) in the collision integrals, determining the kinetic coefficients 
of neutron-star matter [H, [34| . 

Finally, the equations presented in Sec. IV, can be applied to study the response functions and collective modes 
in the equilibrium and in nonequilibrium matter (under the condition th at q vpj <C Aj and u> <C Aj). Furthermore, 
after specifying the collision integral in Eq. (|132l) (see, e.g., Refs. [25l - [30l |54|). the equations of Sec. IV can be used 
to calculate the kinetic coefficients for a superfluid mixture, in particular, the thermal conductivity, shear and bulk 
viscosities. These coefficients are crucial for modeling of the dynamics of neutron stars [5. l6l). 

While doing this work we had in mind possible applications to neutron-star physics. However, the results obtained 
in this paper can be applied to any mixture of strongly interacting superfluid Fermi-liquids, for instance, to ultracold 
Fermi- Fermi mixtures, which have been realized recently [62l [63|. 
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Appendix A 

Let us derive Eq. (|2"3")l for the coefficient on. This coefficient enters the expression (f2"2"| for the matrix A^, describing 
the interaction of quasiparticles with the self-consistent electromagnetic field. Since this matrix is diagonal, it is 
sufficient to consider a mixture of strongly interacting normal Fermi-liquids. The kinetic equation (|19[) then reduces 
to 

- *n» = - ) 6n& + (n£ - n» ) £ P (k,k') 5n^, 

k'cr'j 

+ («i2o-"£o)(«*^-«i^). (149) 
where n^ is given by Eq. © and we use the notation Sn^ = Sn^^ = /a^ ^ a V 
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The last term in the right-hand side of Eq. (|149l) appears due to the interaction of Landau quasiparticles with the 
self-consistent electromagnetic field. The corresponding interaction Hamiltonian has the form 



pal 



kA 



IT I T/ \ (' t (i) 



(150) 



One can easily obtain this term using H em and an equation of motion for the operator ^4 CT aj^ CT (see, e.g., Ref. [24j ] 
In the limit of small q, which is of our primary interest, Eq. (|149[) can be rewritten as 



( \ x (0 On k0 



k'cr'j 



(151) 



The kinetic equation (|151|) [or ()149j) ] is not obviously gauge-invariant. To make the gauge invariance explicit, one can 
notice that k is actually a generalized momentum. It is related to the real momentum p of a quasiparticle i by 



k = p + ej A. 



(152) 



For a very pedagogical discussion of the validity of this expression for Landau quasiparticles, see, e.g., Chapter 3, § 6 
of Ref. [||. 

Thus, we have two momentums, A; and p, and our next step will be to express the distribution function 



pa 



rip g + 5Mpa for quasiparticles with momentum p through the distribution function n k + 5rS ka for quasiparticles with 

(153) 



momentum k. Because of the one-to-one correspondence between k and p, one has 

,(0 -i-JT-W 



J v per 



H k + Sn ko 



In view of Eq. (|152[) . in the linear approximation, 



CO _ (») 

n k0 ~ n p+eiA0 



n po+ Qp e ^ A - 



It follows then from Eq. (|153[) 



a p0 



dp 



a A. 



(154) 



(155) 



Substituting Eq. (|155[) into the kinetic equation (|1 5 1 [) and demanding that the terms which are noninvariant under 
gauge transformations vanish, one obtains the condition 



qvi 



On, 



(i) 
'pO 



Op 



dp 



p'a'j 



III j 



(156) 



or, after performing an integration, the expression (|23[) for oli. For a one-component Fermi- liquid Eq. (|23[) simplihes 
with the help of Eq. ([TJ) and one gets on = e^. 

Using Eq. (|156l) . the kinetic equation (I151[) can be recasted in the well-known gauge- invariant form (see, e.g., Ref. 
45]). 



(uj-qvi) 5n$ = 
where E = i(u;A — qV) is the electric field. 



dn 



(i) 
pO 



<3p 



9 E / y (p,p')*"51'+iei-B 
pVj 



(157) 



Appendix B 

Let us derive Eq. (|64|) directly for a one-component Fermi-liquid. In this case one has to put Jq = and f\ 3 = 
for i ^ j in all equations. In what follows we suppress the particle species indices to simplify notations. We start 
with Eq. (|62[) which can be rewritten as 

A lur =2e^-if l I = 2e^-, (158) 
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where we used the fact that a — e for a one-component Fermi-liquid (see Appendix A) and defined 

I = £*'«n fcV +. (159) 



k' c 



The quantity / is not gauge-invariant. Our aim will be to express / through the vector potential A and the gauge- 
invariant particle current density j, which is given by Eq. (J3HJ) - In view of the definition (|159p . Eq. (|46p can be 
rewritten as 



1 

3 = 2 



m* 2^1 QJ, J 1 k-3, I 



h 3 

hp 



m 



(160) 



Here a and j3 are the space indices. The integral in Eq. (|160[) can be easily taken since the function dNkv/dk has a 
sharp maximum near the Fermi-surface. As a result, we obtain 



1 

3= 2 



1 , fcp/i 



I-e^A. 



m* 37r 2 

As follows from Eq. (JT)) for the effective mass, the expression in square brackets equals 1/m. Thus, one finds 

1 T n A 

J = T^t-e — A 
2m m 

or 

I = 2m (j + e—A) . 

V m J 

Substituting this formula into Eq. (|158[) . one gets 

n ( 1 nfA , a 2m/i , . „ 
A ka = 2e( # )kA - —4± kj = 2e 

\m ki J A: F 



kA c 



Again, using Eq. (J7|) , it follows that the expression in brackets equals 1 / m* . Eq. (|164|) can then be rewritten as 



- kj = 2e 



That is 



2e 2m f 

Aka — — r kA -g 



. . 1 mm* . 

cff = e ~kf fl3 ' 



kA cfi 



(161) 

(162) 
(163) 

(164) 

s 

(165) 
(166) 



This equation coincides with Eq. (|64[) if the latter is written for a one-component Fermi-liquid. 
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